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P ' Abstract 



We investigate scenarios in which a charged, long-hved scalar particle decouples from 
the primordial plasma in the Early Universe. We compute the number density at 
' time of freeze-out considering both the cases of abelian and non-abelian interactions 

and including the effect of Sommerfeld enhancement at low initial velocity. We also 
discuss as extreme case the maximal cross section that fulfils the unitarity bound. 
We then compare these number densities to the exotic nuclei searches for stable 
relics and to the BBN bounds on unstable relics and draw conclusions for the cases 
of a stau or stop NLSP in supersymmetric models with a gravitino or axino LSP. 
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1 Introduction 



The early Universe may have been populated by many exotic particles that, especially if 
charged, should have easily been in thermal equilibrium. No charged relic seems to have 
survived to the present day. In fact there are very strong upper bounds on the density of 
electromagnetically and/or colour charged particles with masses below 10-100 TeV from 
extensive searches for exotic nuclei [1]. The standard lore is therefore that only neutral 
relics may have survived until today. 

However, it is possible that some unstable but very long-lived charged particle froze- 
out from thermal equilibrium and decayed much later to a neutral one. A typical example 
of this kind in supersymmetric models with R-parity conservation is the next-to-lightest 
supersymmetric particle (NLSP) if the LSP and Cold Dark Matter is very weakly interact- 
ing like the axino [2l|3lll] or the gravitino [3, [6]. Recently, such candidates have attracted 
a lot of attention, and indeed the signal of a charged metastable NLSP at colliders would 
be spectacular [UIH]. 

In general, strong bounds on the number density of any metastable relic with lifetime 
of about 1 s or longer are provided by Big Bang Nucleosynthesis (BBN) [0] . They come 
from two classes of processes: on one hand injection of very energetic photons or hadrons 
from decays during or after BBN adds an additional non-thermal component to the plasma 
and can modify the abundances of the light elements fUJ]; on the other hand, if the relic 
particle is electromagnetically charged, bound states with nuclei may arise that strongly 
enhance some of the nuclear rates and allow for catalysed production of e.g. ^Li [U]. The 
bounds of the first type are very tight for lifetimes of the order of 10'^ s and exclude, for 
instance, a neutralino NLSP with a gravitino LSP in the CMSSM An electrically 
charged NLSP like the f can instead escape the first class of constraints in part of the 
parameter space, but it is excluded for long lifetimes by bound state effects [12] • In the 
axino LSP case, the NLSP has a shorter lifetime; the BBN bounds are hence much weaker 
and both, neutralino and stau, NLSP are still allowed [2j. 

In this paper, we investigate the most general case of a scalar charged thermal relic. 
We compute the number density and compare it to the bounds on exotic nuclei for stable 
particles and the BBN constraints for unstable ones. Similar studies have been carried 
out model-independently many years ago [131 [H [IS] for stable relics and we will update 
and improve these computations^ We mostly consider the role of the gauge interaction 
for two main reasons: i) the annihilation into gauge bosons is often the dominant channel 
for a charged particle and ii) it depends only on very few parameters, just the mass of the 
particle and its charge or representation. It is also enhanced by the Sommerfeld effect [T7] , 
analogous to heavy quark production at threshold, which has previously been considered 
for dark matter annihilations in [18l [19], [16], [20l [21] and recently also in the context of 
leptogenesis in [22] . We discuss this Sommerfeld enhancement for the general abelian and 
non-abelian cases. Moreover, we compare the cross sections with the unitarity bound and 
update the unitarity limit on the mass of a stable relic. 

Our main goal is to determine if it is at all possible to evade completely either the 
exotic nuclei bounds or the BBN ones and how strongly the particle has to interact in 
this case. We then apply our findings to the Minimal Supersymmetric Standard Model 
and discuss in more detail the cases of the stau and stop NLSP. 

^Recently the case of general EW charged rehcs as DM was also considered in full detail [IB] , 
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The paper is organised as follows. In Section 2, we briefly review the computation 
of the number density from thermal freeze-out. The formulae for the annihilation cross 
section of a charged particle into gauge bosons are given in Section 3. Here we dis- 
cuss abelian and non-abelian cases, the Sommerfeld enhancement and the unitarity cross 
section. Moreover, we compare the thermal averages with the flrst order in velocity ex- 
pansion. The resulting relic density is discussed in Section 4. In Section 5, we review the 
constraints on stable and unstable relics. These are then applied in Section 6 to the con- 
crete examples of relic staus and stops. Section 7 flnally contains our conclusions. Details 
on the computation of the annihilation cross section and the case of massive gauge bosons 
are given in the Appendices A and B. 



2 Number density of a thermal relic 

The number density of a stable or quasi-stable thermal relic is determined by its annihi- 
lation cross section. In fact the number density of a particle in a thermal bath and an 
expanding Universe is described by the Boltzmann equation 



where the dot indicates the time derivative, C denotes the collision integral of all processes 
that change the particle number and fx is the phase-space density for the particle X. 
For a particle with a conserved parity, like R-parity, the lowest order processes to be 
considered in the collision integral are just two particle scatterings, i.e. annihilations and 
coannihilations. If there is a lighter particle carrying the conserved parity number, C 
includes also the decay into this lighter state, but we will assume that such a decay rate 
is so small it can be neglected at the time of freeze-out and becomes effective only much 
later. Then we have effectively a two step process and we can treat freeze-out and decay 
separately. This is a general feature if the decay takes place via a non-renormalisable 
interaction and is suppressed by an intermediate or even the Planck scale (see e.g. the 
axino [21 H] and gravitino cases [HI E] ) . 

Taking into account only the annihilation of particle and antiparticle, we can write 
the collision integral as [24] 



where (Jann denotes the unpolarised annihilation cross section of an XX pair summed over 
initial and final states. We are here assuming that CP is conserved and no asymmetry 
exists between nx and njf. Note that the production cross section is taken into account by 
the term proportional to fxf^ since we are assuming that the products of the annihilation 
are much lighter than X and are still in thermal equilibrium. 

In this paper we will consider charged relics and concentrate therefore on the annihi- 
lation into gauge bosons, which is the dominant channel in most of parameter space and 
does depend only on the mass and charge of the relic. Note that adding more channels 
only increases the cross section and reduces the relic particle number density further. 
Instead, the inclusion of coannihilations for a charged particle does not always reduce the 
number density as discussed in 
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We can rewrite eq. ([T]) by changing variable to Yx = nx/s, where s(T) = gs^T^ 
is the entropy density, so that the dilution due to the expansion of the universe cancels 
out in the ratio as long as entropy is conserved. It is also convenient to replace the time 
variable with x = thanks to the relation dt = We thus obtain 

- -W.frl-yi) (3) 



15 \gpl mx ^ 



Here we have used H = ^dpjjr, for Mp = 2.43 x 10 GeV, valid during the radiation 
dominated era. Moreover, we define the thermally averaged cross section a^ 



1 f x\ 

where Ki{z) are the modified Bessel functions of order i, characteristic of Maxwell- 
Boltzmann statistics (we are assuming that we can approximate Bose-Einstein statistics 
with Maxwell-Boltzmann statistics). In this expression the rescaled cross section a is 
given by the annihilation cross section averaged over initial and summed over final states 
and multiplied by a factor proportional to the squared M0ller velocity, 

5- = (s-4m^)(T(mx,s) . (6) 

Note that in the centre-of-mass system the M0ller velocity is equal to the relative velocity 
between the annihilating particles and given by 



t;M0i = 2/3 = 2yi-— ^ . (7) 

The rescaled cross section a defined above is dimensionless and function only of x/z = 
^x/v^ (or P) for the case of annihilation into massless gauge bosons and it always 
vanishes at threshold. Then it is easy to see that since we integrate in both x, z, the main 
dependence on the charged relic mass is contained in the prefactor in eq. (jll) and can 
be reabsorbed in a rescaling of Yx Yx/mx- For this reason we obtain nearly exactly 
Yx oc mx if there is no other mass scale involved. Note that in principle a much weaker 
logarithmic dependence on mx is present in the value of the freeze-out temperature, when 
Y begins to deviate from Y^g. 

We are here computing the yield of the particle X and to obtain the yield of particle 
and antiparticle we multiply by a factor of 2 or divide the cross section by 1/2, since we 
are assuming nx = Ux- Also note that, contrary to intuition, for a particle with internal 
degrees of freedom like a coloured state, the total yield is the solution of the Boltzmann 
equation (jl]) with the cross section averaged over the initial states. Instead the yield per 
degree of freedom is obtained from the cross section averaged over X, but summed over 



^Note that our definition diff'ers from the one in by a factor m\/x^ since we prefer to work with 
a dimensionless quantity and to absorb here aU the dependence on x. 
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The presence of many degrees of freedom in the initial state has then the effect of 
partially compensating the large cross section coming from the multiplicity of the final 
states. 



3 Annihilation cross section for a charged particle 
into gauge bosons 

3.1 Abelian case 

For an abelian gauge symmetry, there are only three Feynman diagrams contributing to 
the annihilation cross section, analogous to those shown in Fig. [1) the t- and u-channel 
exchange of the scalar particle itself, and the 4-boson vertex. The amplitude is symmetric 
in the exchange of the gauge bosons and for a particle of charge exgi it is given by 
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The cross section is a function of the mass and charge of the relic: 
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(9) 



where ai = gf/^iir) is the gauge coupling; note that a symmetry factor 1/2 has to be 
added due to the symmetric final state of identical particles. For the rescaled cross section 
this gives 



2^ ^ 4/? 
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log 
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1 + /?. 



which is a function only of /? = a/1 — Amj^/s and the charge of the particle. 



(10) 



3.2 Non-abelian case 



The computation for the annihilation into non-abelian gauge bosons is slightly more 
involved, since there is an additional contribution from the Feynman diagram with a 
gauge boson in the s-channel and the 3-gauge-boson vertex. The amplitude can be divided 
into a symmetric and an antisymmetric piece in the group indices. The symmetric one is 
analogous to the abelian case: 
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In fact any rescaling of the cross section by a factor p due to a different counting of the degrees of 
freedom can be absorbed into a rescaUng 1/p of the yield(s). 
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Figure 1: Feynman diagrams for the annihilation into gauge bosons, here for the case 
of gluons. In the abehan case, there is no 3-gauge-boson vertex, so the last diagram is 
absent. 



while the antisymmetric part is given by 
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(12) 



The two contributions do not interfere due to the different symmetry, so we have for the 
amplitude squared, summing only over physical polarisations of the final gauge bosons: 
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(13) 



Then the sum over all final and initial states for a scalar in the fundamental repre- 
sentation T° of the gauge group SU{N), normalised such that Tr(T"T^) = 5"'^ /2, can be 
obtained from the usual group invariants: 



j,i,a,b 



X(iv)(i.l(cliiv)-4))^ (^^-^)f^ 



(14) 



where we have separated the singlet and adjoint contributions to the symmetric part 
for later convenience, included a factor 1/2 for identical particles in the final states and 
used the Casimir invariants for the fundamental and adjoint representations, Cf{N) = 
, Ca{N) = N. Note that the ratio of the singlet to adjoint contributions is given 
simply by tjtzi- The antisymmetric channel instead gives 
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(15) 



Finally we obtain for the cross section averaged over initial states: 
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(16) 



This result coincides for = 3 with that reported in [26j for the Born cross section of a 
pair of gluons into squarks, allowing for the exchange of the initial and final state. 
Then the rescaled cross section for SU{N) is 



2'Ka 




(17) 



Note that the contribution of order (3 in the expression above in the limit /5 — > is due 
to the symmetric part of the matrix element and that the antisymmetric piece instead 
vanishes at that order. Therefore the symmetric part of the cross section dominates at 
threshold. 

So we see that for a non-abelian interaction the cross section is larger than for the 
abelian case, not only due to the possibly larger coupling a^r, but also due to the opening 
of an antisymmetric channel and of course to the multiplicity of the final states. In fact 
for large N the averaged cross section increases as N and therefore the yield decreases as 
l/N. 



3.3 Annihilation into SU{N) gauge boson and photon 

The annihilation cross section into gluon and photon is just the same as the abelian one, 
but with a different vertex for the gluon. Then considering a particle of electromag- 
netic charge exgi, in the representation of the gauge group SU{N) with coupling g^, 
annihilating with its own antiparticle, the amplitude is given by @ 
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From this we easily obtain the cross section as: 



(TiN{mx, s) 



t: 



« |2 
ji\ 



2m\' 



(i- 



log 



v + 



1 + 



Am?x 



1 _ l!!^ 



4m, 



(19) 



Strictly speaking, in this case the final state particles are diff'erent and therefore there are no inde- 
pendent t- and u-channels, but we can still write the amplitude to be symmetric in t and u in order to 
make direct contact with the previous results. 
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where ai^Ar are the gauge couphngs and the symmetry factor 1/2 in this case is absent 
since the final particles are not identical. 

Averaging over the initial and summing over the final state, we have 



N 



2N 



for the fundamental representation. This gives for the rescaled cross section 
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which is a factor (A^^ — l)aN /{Naiej^) larger than the pure U{1) contribution. Again the 
cross section increases as A^ for large A^. 



3.4 Annihilation into physical Z and SU (N) gauge boson/photon 

The annihilation cross section into massive Z and photon/SU(N) gauge boson has the 
same form as the abelian one. We consider here a particle with Z-coupling giCz, in the 
representation T" of the gauge group SU{N) with coupling g^, annihilating with its own 
antiparticle and we obtain 
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[22] 



where p4 is the Z boson momentum obeying p\ = M|; the annihilation into photon and 

g[ex- Then we easily obtain the cross section 
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as: 
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(23) 



where ai^N are the gauge couplings. 

Averaging over the initial and summing over the final state as in eq. ([20 
rescaled cross section 



we have for the 
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where az = Ml/m?x. Note that the cross section for annihilation into photon and Z, 



is given by the substitution 



2N 



a[ex- For the specific case of the right-handed 



stau (stop), the coupling with the Z boson and photon are respectively given by e|ai 



aem^SiVr 6w (e|ai = A/Q af^^iaxr 9\f^r) and exd'i 
the Weinberg angle. 



A/9aem), where 9w is 
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3.5 Annihilation into massless EW gauge bosons 



The cross section for annihilation into massless SU{2)l gauge bosons can be obtained 
directly from the general formula for the non-abelian case. One has to take into account, 
however, that in this case the scalar SU{2) l doublet is not degenerate in mass and that 
the initial particles can be a mixture of left- and right-chiral states. We neglect here the 
effects of EW symmetry breaking; the results are hence applicable for the case of a heavy 
relic that decouples before EW symmetry breaking takes place. 

Considering the scalar relic to be X = Xl cos 6 + X/j sin 6 and denoting with mx' the 
mass of its left-handed doublet partner, which is sufficiently larger than mx to neglect 
coannihilations, we obtain for the annihilation cross section into W^'"^ gauge bosons: 

aw2is,mx,mx') = 
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s — 4m^ 
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(25) 



while the annihilation into is similar to the abelian one in eq. ^ for ex = cosd/2. 
Note that the cross section is suppressed by the mixing angle as cos^ 9 and by the fact 
that the group indices are not summed for the initial state. Also in this case the rescaled 
cross section is not just a simple function of /?, but also of the mass difference in the 
doublet. We have in fact 
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(26) 



where 5"^ = {mj^, — m\)/s. The cross section still vanishes for /3 = and is finite for 
52 — i> 00. The detailed expressions for the case of broken EW symmetry are much more 
involved and include also the contribution of the Higgs s-channel allowing for resonance 
enhancement. They are given in Appendix B. 



3.6 Sommerfeld enhancement 

In the previous sections we have computed the annihilation cross sections to lowest order 
in the gauge coupling. However, it was shown long ago [T7j that an expansion in terms 
of the coupling is inadequate close to threshold, where the velocities of the annihilating 
particles go to zero, 

/J.fTM^O. (27) 

The enhancement at low velocities becomes apparent when one computes the one-loop 
corrections, which are enhanced by a factor Here, C is a process-dependent constant, 
a is the gauge couphng of the annihilating scalars, ai in the case of U{\) boson exchanges, 
or a^v for SU{N) gauge boson exchanges, respectively. To account for this long-distance 
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effect, one therefore has to resum a whole class of diagrams, which consist of t-channel 
ladder-type exchanges of massless soft Coulomb SU{N) or U{1) gauge bosons between 
the annihilating charged particles. 

This resummation of terms ~ a^/jS^ leads to the so-called Sommerfeld factor which 
multiplies the lowest-order annihilation cross section. The Sommerfeld enhancement is 
given by the modulus squared of the particle wave function at the origin, 

E^mO)f= , \ y ^ = (28) 

1 — exp(— [j 

Because this effect is a long-distance one, taking place at a scale ~ Pmx, it factorises 
from the annihilation cross section which is a short-distance effect at the hard-scattering 
scale of order of the mass rnx- Schematically, 

a'^iP, mx) = E{a{Pmx)) x . (29) 

Here, cr'' is the leading-order annihilation cross section, which has been presented in the 
preceding subsections. Eq. fl^Ul) is in principle only valid if the annihilating partons are 
in a single SU{N) channel, i.e. for particle in the fundamental representation either in 
the singlet or adjoint configurations. If multiple channels c contribute, eq. (l29l) has to be 
modified to 

iP,mx) = Y^EcWmx)) x . (30) 



Here, cr°(/5) is the projection of the leading-order annihilation cross section in the rele- 
vant channel. For a scalar in the fundamental representation of the SU{N) gauge group 
annihilating into massless SU{N) gauge bosons, we have seen that only the contribution 
proportional to the group-symmetric part survives in the limit of vanishing (3 and is en- 
hanced at low velocities. Therefore at leading order the cross sections can be taken to 
be the same for the singlet and adjoint part up to colour factors and proportional to the 
total cross section given in eq. (ITTll 0. We note also that due to the presence of more than 
one channel, the Sommerfeld factor for an SU (N) gauge theory becomes dependent also 
on the final states, since not all channels may contribute to the annihilation into a given 
final state. 

However, the presence of the thermal bath complicates things, as the interactions with 
the background gauge bosons may prevent the annihilating partons to be initially in a 
definite SU{N) channel. The time scales for the Sommerfeld effect and the interactions 
with the thermal bath are of competing order, so it is not clear how strong such effect can 
be. In this paper we will consider both extreme situations, i.e. the case when the thermal 
bath has no effect and the case when there is no definite initial channel. In the latter 
case, it was argued in the literature that due to the mixing of states one should just take 
an average C"*" extracted from the averaged one-loop correction, leading again to a single 
Sommerfeld factor as in eq. fl2^ (see, for example, ref. [2Z!)- While the two approaches 
give identical results by construction at first order, they correspond to two quite distinct 
resummations of the higher orders and they are numerically substantially different. 

We obtained the coefficients C by computing the l//3-enhanced contributions for t- 
channel SU{N) gauge boson exchange at one loop in the threshold expansion (see for 



^ Taking the true crj a.nd cr^ instead, differs from the total ct'^ only in the terms suppressed by and 
amounts to a correction smaller than 1% at threshold where the Sommerfeld factor is effective. 
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example [2H] and references therein). For the generation of the relevant one-loop graphs 
and the Lorentz algebra we used the Mathematica packages FeynArts and FeynCalc [29] . 
We simplified the resulting expressions to only keep terms that are leading in j3, that is, 
we only kept terms that are enhanced in the soft region of the one-loop integrals, which 
were then simple enough to perform by hand. Alternatively, as mentioned above, one 
obtains the form (^E\i directly by computing the normalised wave function at the origin 
from the Schrodinger equation, describing the annihilating parton pair, with a Coulomb 
interaction potential for positive energies ~ P'^mx 

The Sommerfeld enhancement due to exchanges of massless Coulomb SU{N) gauge 
bosons is the same for the singlet channel of annihilation into SU{N) gauge bosons Bi^ 
and the annihilation into U{1) gauge bosons Bi , 

The factor for the adjoint channel is instead found to be negative and thus suppressing. 



Cts^B^B, = CriN) - = . (32) 



The same factors or apply also for other final states of the singlet or adjoint 
channels. For example, the Sommerfeld factor for tt* hh is Cg^(^^■^ = 4/3, while that 

for it* gh,g-f,gZ is C^u{3) = -1/6- 

Even if the adjoint channel leads to a suppression, upon summing over both contri- 
butions in eq. (l30ll . the net effect is still quite enhancing for small N. We have then in 
fact 



E^{a{pmx)) X — — - + EA{a{pmx)) x 



2 



(33) 



where, as described above, we have taken (iV^ -2)/2 = {N'^ - 2) / {N'^ - 4) cr^ = (t°(/5), 
and cr^{P) is given in eq. (ITTll . For SU{3) this gives 



4miP.nix) = 4uis)i/3)^ 



16 



4 TTQ^ 1 TTQtg 

e~3 13 1 _ e6 /3 



(34) 



so that the enhancement in the singlet dominates over the suppression in the adjoint 
channel. 

On the other hand, averaging the one loop contribution over initial channels results 
in a factor 

A^^ + 2 
2A^(A^2 _ 2) 

which is although enhancing, much less so than the net effect of the summation over 
singlet and adjoint channels. For SU{3), this factor is CH^^g-, = 11/42, leading to 



Css^BnBn - ^7^77^72^79^ ~' ^suiN) ^ (35) 



<ui3)W, mx) = a°^(3)(/3) ^ (36) 

1 — e ''^/s 



^Averaging over initial channels is not to be confused with averaging over initial states which is to be 
done in addition when solving the Boltzmann equation. 
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— m = 100 GeV, SU(3) 

— m = 1 TeV, SU(3) 

— m = 500 GeV, SU(10) 



"O.OO 0.05 0.10 0.15 0.20 0.25 0.30 

Figure 2: Ratio of summed over averaged Sommerfeld enhancement, <y^su{^ / '^TuJ.?,)-! ^ 
function of (3. The full red line shows the SU{?)) case for a mass m = 100 GeV and the 
dashed blue line for a mass m = 1 TeV; the dotted green line is for the hypothetical case 
of SU{10) with m = 500 GeV. 

Note that the first term of the expansion of eq. fl36l) coincides with the 1-loop result of 
[26] for gg —>■ qq* near threshold. 

If the difference in the exponents in the denominators of eqs. (1341) and (l36l) could be 
neglected the two expression would be equal. However, in the small (3 region where the 
Sommerfeld enhancement is relevant, the difference amounts to up to 50 % for 5'?7(3), 
and is even larger for hypothetical larger iV, see Figure |21 

For scalars charged under a U{1) group, there is of course a corresponding enhancement 
due to U{1) boson exchanges. However, the enhancement factor is now governed by the 
U(l) coupling, and thus weaker than an enhancement under a strong SU (N) gauge group. 
The Sommerfeld factor for the dominant annihilation channel into U{1) gauge boson pairs 
can very simply be determined from the abelian part of the calculation that led to the 
factor quoted above. We find, 

CsS^BiBi = 1 ! (37) 

for t-channel U{1) exchange, and the coupling in eq. (l28i) is the U{1) coupling ai. 

Another issue regarding the thermal bath is the fact that gauge bosons acquire a mass 
through interactions with the plasma. This Debye screening effect happens at a scale 
of order ~ gT, whereas the Sommerfeld effect is of order ~ amx(3 ~ ay/rrixT ^ gT. 
Thus the thermal masses of initially massless gauge bosons do not affect the Sommerfeld 
enhancement. 

Finally, there are also massive gauge bosons such as Ws and Zs to consider. The 
Sommerfeld factor arises from instantaneous Coulomb exchanges of massless gauge bosons 
between the slow moving annihilating pair close to threshold, thus resulting in an 1//3 
enhancement, signalling the inadequacy of trying to describe this exchange in an expansion 
in terms of loop corrections. Naturally, massive gauge bosons have a finite width, and 
thus cannot be exchanged instantaneously. In terms of Feynman graphs, the momentum 
flowing through a massive gauge boson that is exchanged between the annihilating pair 
is naturally cut off by the mass of the exchanged boson and can never become too soft. 
The Sommerfeld effect is exponentially suppressed with the mass of the gauge boson, as 
an analysis of the wavefunction picture reveals. It can nevertheless become important for 
relics with masses much larger than the electroweak scale, as a very heavy Wino discussed 
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in [in] . In the following we will consider only the case of massless gauge bosons, which is 
the dominant effect for coloured relics and for purely right-handed sleptons. For a more 
detailed discussion in case of massive EW gauge bosons we refer the reader to [19], [16] . 



^non—el,J , ^ \2 ^ V 



3.7 Unitarity bound 

We next compare the above cross sections with the unitarity bound. Using unitarity and 
partial wave expansion, the non-elastic cross section for a particle with spin Sp is given 
by [E] 

47r (2J+l)(l- r^5) 
{2sp + 1)2 

where J is the angular momentum of the process, Pi is the initial particle momentum, 
4p? = s/?2 in the centre of mass frame in our case, and rjj is the contribution of the 
elastic part. This gives an upper bound for the annihilation cross section with angular 
momentum J as 

167r(2J + l) 

The lowest value is obtained taking J = and since the s-wave annihilation is usually the 
dominant contribution for a scalar non-relativistic particle with Sp = 0, we will take it as 
a reference value. We therefore have for the maximal rescaled cross section: 

^max = lQn (40) 

independent of the particle mass or energy. In this case the thermal averaging is simple 
and we obtain 

167rfs:2(2a:) 

which we will consider in the following to be the maximal cross section per degree of 
freedom^. We see clearly that the cross sections discussed above satisfy this bound and 
are suppressed at the very least by a^. Figure [3] shows the rescaled cross sections for 
the abehan and non-abehan cases, eqs. (fTUj) and (|T7I) . together with the unitarity bound 
eq. fl40l) as a function of the relative velocity of the annihilating particles. 

The unitarity cross section dmax can be used to obtain a lower bound of the yield. 
Moreover, it can be taken as the maximal annihilation cross section possible even after the 
QCD phase transition, when the coloured states are confined into the equivalent of scalar 
hadrons and fermionic mesons [30] . Constraints from cosmology on such kind of hadronic 
states have been mostly studied for the case of a stable exotic quark [H], a gluino LSP [27] 
or for very long-lived gluino in the split SUSY scenarios [18]. It has been argued in [31] 
that the annihilation cross section for such states could become much stronger, if bound 
states between two scalar hadrons/fermionic mesons are formed with rate ~ vt/Aqj^^, 
and in that case the coloured relic abundance after the QCD phase transition is further 
reduced below Y ~ 10^^^ — 10^^''. We will not consider this possibility in the following, 
but note however that, while most of the cosmological bounds for a decaying relicare then 
satisfied, one still needs to consider the bounds for a stable relic. 



\Pmax'V)x = / N2 ; (41) 



Note that here we are computing expHcitly in the centre of mass frame, while the Bohzmann equation 
requires to use the covariant or lab frame. The difference between the two frames has been discussed in 
|24j and gives only a small correction for non-relativistic particles, which we neglect here. 
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Figure 3: Dependence of the rescaled cross sections on the relative velocity /?, nor- 
malised to 1 at large s, i.e. f3 = 1. The solid lines show the leading order section, 
the dashed/dashed dotted lines the effect of the Sommerfeld enhancement, that makes 
the cross sections non-vanishing at the threshold (3 = 0. The SU{3) cross sections are the 
upper (red) lines, including the averaged Sommerfeld factor in the dashed line and the 
summed one in the dash-dotted. For the abelian case (blue lines) the Sommerfeld effect is 
much milder and shown in the dash-dotted line. Note that the region for /3 ~ contributes 
more strongly to the thermally averaged cross section due to the Boltzmann-suppression 
for large j3. 



3.8 Thermally averaged cross sections and velocity expansion 

We integrate eq. ([5]) numerically to obtain the thermally averaged cross section. Very often 
such a quantity is instead approximated with the first terms of its velocity expansion, since 
the relevant regime takes place when the annihilating particles are already non-relativistic. 
To obtain such an expansion, one can use the approximation 

s - Amx ^ 4m^/52 (42) 

and expand in /5 the expression 

a.M.i-^a(/5) . (43) 

We see that if a is constant at zero velocity, the cross section is enhanced like 1/(3 in that 
limit. This is indeed the case both for the Sommerfeld-enhanced cross section and the 
unitarity one. 

The first term in the expansion, which is independent of the velocity and coincides 
therefore with the first term in the expansion of the thermally averaged cross section [21] , 
is given by 

aa,v - ^^^ + 0{(3'), (44) 



m 



X 



7ra% (iV^-l)(iV 2 - 2) 
ml Am 



^na,v '- + 0{P^), (45) 
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Figure 4: Ratio of the thermally-averaged cross section and the first term in the velocity 
expansion around /? = 0, for mx = 350 GeV. The thick solid line is for the abelian, 
the thin line for the non-abelian {SU{3)) case. Dash-dotted and dashed respectively are 
the same ratios including the Sommerfeld enhancement, only the averaged one for the 
non-abelian case: we see that in this case the thermally averaged cross sections do not 
converge to the first order term in velocity, but that the latter can still give a good estimate 
within 15% of the full result in the abelian case; for the non-abelian case the Sommerfeld 
enhancement changes the result considerably and the velocity expansion fails. Note that 
the case of the summed Sommerfeld factor is outside the range of the plot. 

for the abelian and non-abelian cases respectively. 

We plot in Figure Hlthe thermally averaged cross sections as a function of x normalised 
with respect to the first term in their velocity expansion including also the Sommerfeld 
enhancement factor, both for the abelian case and for the QCD case with = 3. We 
see that keeping only the lowest order overestimates the thermally averaged cross section, 
i.e. underestimates the yield, in the abelian case by at most 20% in the region of freeze- 
out {x ~ 30). The non-abelian case for = 3 is approximated better also because 
the freeze-out takes place at a larger x ~ 40, i.e. smaller f3. On the other hand, once 
we include the Sommerfeld enhancement, the thermally averaged cross section does no 
more converge to the first constant term in the velocity expansion due to the threshold 
singularity at /? = 0. Nevertheless the first order term without the enhancement can 
still give a reasonably good approximation for the abelian case, since the Sommerfeld 
enhancement partially compensate the 20% underestimation of the Born result. For the 
non-abelian case the Sommerfeld enhancement is so strong that the low energy expansion 
can give only an order of magnitude estimate. 

4 Results for the relic density 

We solve the Boltzmann equation (jlj) numerically for the exact thermally averaged cross 
sections given above. This improves the old results ^13j that were obtained with the 
velocity expansion. 

For the case of an abelian charged relic, we consider ex = ±1 and we set the coupling 
to be ttem = 1/128. For the non-abelian case we take = 3 and oa? to be the QCD 
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Figure 5: Time evolution of the particle yield for the cases of abelian and non-abelian 
cross section, for mx = 200 GeV. The upper (blue) curve is for an electromagnetically 
charged scalar particle with unit charge, while the lower (red) curves correspond to a 
single coloured scalar in the fundamental representation without the Sommerfeld factor 
(solid) and with the Sommerfeld factor averaged (dashed). We see that the treatment of 
the Sommerfeld factor has an impact of about 30% on the final number density. 




coupling (Q) with Q = 2mx in the hard process and Q = (3mx in the Sommerfeld 
correction, c.f. Sect. 13. 6[ In order to avoid the non-perturbative regime, we cut off the 
running of as at Q = 2 GeV, i.e. a3(Q < 2 GeV) = as (2 GeV). 

For the entropy and energy density parameters we take g^J"^ = = 10, since we 
expect the freeze-out to take place between 10-100 GeV, when only the light Standard 
Model particles are still in equilibrium in the thermal bath. 

Our results are plotted in FigureO We see that the yield Y follows relatively closely the 
equilibrium density until the time of freeze-out, which happens at different values of x for 
the different cross sections. As expected the non-abelian interactions being stronger gives 
a considerably lower relic density. The ratio between the two cases is well approximated 
by the ratio of cross sections, (Tnabv/cTabV, at zero velocity : 



Yah _ 7 aj 
Y„,„, 27 a2 



40 . (46) 



em 



We next consider the dependence on the only dimensional parameter, the mass of 
the charged relic. We have seen that the thermal average can be written only as a 
function of x and since we are integrating the Boltzmann equation to x — *• oo we get 
rid of the dependence on mx that is contained there. A subleading dependence would 
survive by integrating to a finite value of x, but this effect is negligible for the present 
universe with a temperature Tnow ~ lO^^eV <^ mx- On the other hand, the mass 
directly enters in the coefficient of eq. (jl]) and that is the stronger dependence on mx- 
Note that this dependence is present even in the unitarity case, where the reduced cross 
section is explicitly independent of the mass and velocity. In general therefore the yield 
is proportional to the mass and can be rescaled as 

y(mx)=y(lTeV)(^) . (47) 
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Figure 6: Dependence of the yield on the mass of the charged rehc. From top to bottom, 
the first (blue) line is for the case of an electromagnetically charged relic, while the second 
(red) line is for a coloured relic, the dashed and dash-dotted lines include the Sommerfeld 
factor, averaged and summed respectively. The lower two (black) lines correspond to the 
maximal annihilation cross section given by unitarity - the solid one for a single d.o.f., 
the dotted one for 3 d.o.f. for the fundamental representation of QCD. Note that the 
non-abelian case is still three orders of magnitudes away from the unitarity cross section. 



with F„f,(lTeV) = 3.9 x lO^^^ ^nd r„afc(lTeV) = 1.6 x lO^^^ for the abelian and non- 
abelian cases, respectively, for the total degrees of freedom, including antiparticles. For the 
case of the unitarity cross section, the total yield becomes instead Yumil TeV) = 6.6x 10~^^ 
(or 2 X 10"^'' for three degrees of freedom). 

Since the energy density also increases for larger masses, this can be used to give a 
constraint on the mass of any stable thermal relic from the maximal cross section allowed 
by unitarity [15]. Using the WMAP 5- year results [32] for the most conservative upper 
bound for the matter density, we can update such bound. In fact imposing 



nxh"^ = mxYx+x{Tnow)s{Tnow)/ Pc < 0.13 
gives us for a single degree of freedom the constraint 



mxYx+x{Tn 



GeV 



< 4.6 X 10 



-10 



resulting for a scalar particle in 



mx < 280 TeV . 



(4J 



(49) 



(50) 



Note that for a fermionic spin 1/2 relic the unitarity cross section is reduced by a factor 
four and therefore the bound on the mass is stronger by a factor two. 



5 Constraints on cosmological relics 

We review here the constraints on the abundance of cosmological relics that we will 
compare with the number density of a charged scalar relic in the next section. First we 
will consider the case of stable relics (i.e. with lifetimes longer than 10^^ s) and next relics 



18 



with lifetimes in the window 0.1 — lO^^s. Note that for shorter hfetimes the constraints 
are non-existent, as long as the particle did not dominate the universe dynamics before 
decaying or produce a large amount of entropy, while for lifetimes between 10^° — lO^^s 
bounds from CMB distortion [33] and from the measured photon diffuse flux [3l] apply, 
but will not be discussed here. 



5.1 Stable relics 

The possibility of existence of some more exotic cosmological relics than the known light 
elements stimulated many years ago the search for exotic nuclei in water and other ma- 
terials on the earth. Those searches were unsuccessful and provide a very strong limit on 
the number density of any relic that would bind electromagnetically with an electron or 
in nuclei, under the assumption that such particles are equally distributed in the Universe 
compared to baryons. If such relics were present long before structure formation, it is 
highly probable that they were trapped together with baryons when the universe's density 
was still nearly homogeneous, so that we can expect their number density not to be too 
strongly dependent on the local environment. Note that in any case these bounds are so 
strong that the possibility of such a relic to be Dark Matter is completely excluded. 

The most recent constraints are those obtained by [35] looking for anomalously heavy 
hydrogen in deep sea water, which apply to an electrically positively charged relic, and 
give for masses 5 GeV < rrix < 1.6 TeV: 

Fx+ < 4 X 10-1^ Yb = 3.5 X IQ-^' ( (51) 

Taking into account the gravitational effect in deep sea, this corresponds to a concentration 
of the order of 10~^^ at sea level or equivalently 

Yx+ < 0.9 X 10-38 ( \ , (52) 

which is comparable to other limits in the same mass range, [36]. For larger masses up to 
a TeV, a slightly looser bound Yx+ /Yb < 3 x 10"^° was found by [37], while for even larger 
masses 10 TeV < rrix < 6 x lO'* TeV it weakens even further to Yx+/Yb < 7 x 10~^^, as 
given by [SB], i.e. 

rx+ < 6 X 10"^^ ( . (53) 
V0.0223/ ^ ^ 

For electromagnetically neutral, but coloured relics, the bounds are obtained from 
considering heavier elements and are considerably weaker; using the results of [3Z] for 
Carbon, the limits are of the order Yx+/Yb < 4 — 8 x 10~^° for mx = 0.1 — 1 TeV, reaching 
2 X 10~^^ at the largest mass considered 10 TeV. For larger masses rrix < 100 TeV only 
the constraint by for lead is present, giving 

Fx < 1.5 X 10-13 Yb = 1.3 x 10-^^ ( . (54) 

V0.0223 / ^ ^ 

We see that these constraints are very strong. In order to reach even the weakest bound 
of Yx < 10~^^, the unitarity cross section is way too weak and needs to be increased at 
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least by nine orders of magnitude, i.e. 

5](2J+l)>10^ (55) 
J 

Therefore stable relics are allowed only if their interaction does not belong to the Standard 
Model and they cannot form exotic atoms/nuclei or if their annihilation rate becomes 
much larger than the unitarity one as it can happen if they interact strongly and can 
form intermediate bound states. But in any case, note that cross sections of the order 
tt/Aq^^ that can arise after the QCD phase transition are not sufficient to evade these 
constraints [31], so their interaction would have to be stronger than QCD. 

5.2 Unstable relics 

Different cosmological constraints exist on the density of an unstable relic, depending on 
its lifetime. For lifetimes between 0.1s and 10^° s, the strongest constraints come from 
Big Bang Nucleosynthesis. In fact, if the relic decay injects very energetic particles into 
the thermal bath during BBN, it can change the abundances of the light elements. Since 
standard BBN agrees quite well with the primordial abundances of Helium-4, Deuterium 
and (within a factor of two) Lithium-7 inferred from present astronomical observations , 
the relic density has to be low enough not to change those predictions too strongly. These 
effects are present for any decaying particle and have been studied in various papers 
(see p!Ol HOl HU |12] and references therein). For lifetimes above 3000 s, corresponding 
to the time of production of Lithium, additional constraints are present if the relic is 
electromagnetically charged and can form a bound state with positively charged nuclei 
increasing the rates for Lithium-6 production [TT]. The Standard BBN prediction for 
the ^Lz abundance is actually way too small compared to the observed one, so that the 
presence of a charged relic with appropriate lifetime can help reconciling BBN with the 
measured abundances of ^Lz, ^Li [13], but we will disregard this possibility and only 
concentrate on the exclusion region. 

We summarise here the main results from various BBN analyses and give conservative 
bounds on the energy density of the decaying relic and compare them with our com- 
putation of the relic density. Since we are interested in escaping the BBN constraints, 
we focus mainly on the strongest bounds, but we keep conservative values for the light 
element abundances. Note that in many of the analysis slightly different ranges for these 
abundances are considered, corresponding to shghtly different constraints on the decaying 
relic. 

In general, the decay can produce very energetic SM particles that can initiate either 
hadronic or electromagnetic showers in the plasma. The most stringent bounds are ob- 
tained for a relic that produces mostly hadronic showers, since electromagnetic particles 
like photons or electrons can thermalise very quickly by interacting with the tail of the 
CMB distribution until times of about 10^ s. So we will consider in the following the 
constraints for relics producing hadronic showers with a branching ratio Bh = 1- We 
will comment later on the case where this branching ratio is smaller. There are then 
practically three regions of the lifetimes as discussed in |40j : 

• 10^^ s < r < 10^ s : the dominant effect is the interconversion between protons and 
neutrons, that changes the Helium abundance, overproducing it; 
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• 10^ s < T < 10^ s: hadro dissociation is the most efficient process and the bound 
come from the non-thermal production of Li and D\ 

• 10^ s < r < 10^^ s: photodissociation caused both by direct electromagnetic show- 
ers and by those generated by the daughter hadrons starts to dominate and the 
overproduction of ^He is the main result. 

It is clear that these limits depend on the decay branching ratio Bh into hadrons for 
lifetimes r < 10^ s, while they are independent of Bh for longer lifetimes. In Table [1], 
we give conservative bounds taken from the general analysis of [40J for the three regions, 
assuming B^ = 1. Similar constraints were obtained independently also by [41j. Note 
that the bound for short lifetimes becomes approximately one order of magnitude weaker 
if one takes a more recent value of the ^He abundance as discussed in [32j. Unfortunately 
this new publication does not provide constraints for a general relic, but discusses only 
the explicit cases of a bino neutralino or a right-handed stau. 
The limits we use can be parameterised as 

Yx+x < 1-0 X 10-'' {i^) ° ' ~ - s , (56) 

Yx+x < 1.1 X 10-16 (^_^^ fo^ ^ _ g , (57) 

The assumption Bh = 1 is surely valid if the decaying relic is coloured, while Bh can be 
different if it is only electromagnetically charged, as in the case of the stau. If the branch- 
ing ratio into hadronic modes for the relic is less than one, the hadronic BBN bounds are 
relaxed accordingly by a factor 1/Bh- For intermediate lifetimes, then electromagnetic 
showers can become a more important effect, but only if Bh < 0.01. 

For electromagnetically charged relics with lifetimes longer than about 3000 s and low 
Bh < 0.1 — 0.01, strong bounds also come from considering the catalysed overproduction 
of ^Li [TT]. In fact when bound states between nuclei and the relic can form such as 
^HeX~ , many nuclear rate are modified and change the final abundance especially of ^Li 
and "^Li. For particles decaying after 5 x lO^s it has been argued that uncertainties in the 
nuclear rates make such constraints weaker than the general ones discussed above so 
we will consider here catalysed BBN constraints only for the intermediate lifetime range. 

Unfortunately, different values for these bounds are given in the literature; in [13, US] 
they are found to be maximally at the level of Yx- < 1.4-2 x 10"^^, while the latest value 
in [44j is maximally Yx- < 10~i^, taking a larger window for the ratio ^Lif^Li. Here we 
will use as a constraint the simple interpolation for the total yield |^ 

Y- < / 2 X 10-^ (3^) ' for rx < lO^s ^^g^ 
\ 2 X 10-15 for Tx > lO^s ^ ^ 

that lies somewhat in between. The bounds from catalysed BBN do not apply for coloured 
scalar relics because these should have a large branching ratio into hadrons, such that the 
'conventional' BBN bounds from hadronic showers are much stronger. In passing note 
also that up-type squarks would mostly hadronise into neutral fermionic mesons which 
are lighter than the charged ones 



^The catalysed BBN constraints restrict only the abundance of the negatively charged particles, but 
we give here the constraint for the total yield assuming 2Yx- ~ ^x+x- 
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Maximal values of mxYx+x (GeV) allowed by BBN 



mx (TeV) 


10-^ - 10^ s 


10^ - 10^ s 


10^ - 1012 s 


0.1 


2 X 10-" 


5 X 10-1^ 


10-14 


1 


1 X 10-1° 


10-13 


10-14 


10 


5 X 10-1° 


3 X 10-13 


10-14 



Table 1: Maximal allowed values of mxYx+x i^i the different region of lifetimes taken 
from Figures 38-40 of [50]. We are assuming here that the energy released in Standard 
Model particles is one half of mx as happens in a two body decay of the NLSP into LSP 
and the NLSP non-supersymmetric partner and that all the energy is released in hadrons. 
In general the strongest bound is for longer lifetimes and it is independent of mx and the 
hadronic branching ratio. The bounds in the second column come from D, but the ^Li 
ones, that are sometimes considered too strong [13], are not very far away. 

We summarise the constraints in Fig. [TJ which shows our conservative bounds in the 
plane of total number density vs lifetime. Note that the constraint from catalysed BBN 
are for the stau stronger than the hadronic ones for lifetimes longer than ~ 10"! s and 
exclude a light stau NLSP with a 100 GeV gravitino LSP in the CMSSM [I2]. 

Comparing with Fig. [6l we see that even for a charged relic that can annihilate ef- 
ficiently, the BBN bounds are very strong; in particular the case of a simple abelian 
interaction seems to be excluded for any charged relic whose lifetime is longer than 0.1 s 
and produces hadronic showers with Bh = 1- For the coloured case the situation is less 
severe, but even with the Sommerfeld enhancement, which reduces the yield substantially, 
it is not possible to evade the bounds completely. Still all masses above approximately 
50 GeV are excluded for lifetimes longer than 100 s, while for shorter lifetimes masses up 
to 700 GeV are allowed. A much larger number of colours than three would be needed to 
relax all bounds. Even the unitarity case reaches the strongest BBN constraint at masses 
around 700 GeV for 3 degrees of freedom or 1 TeV for a single one. 

6 Application to the MSSM 

Until now we have considered the ideal case that the relic particle has only one single 
interaction. In realistic models, however, more than one interaction - and hence more 
than one annihilation channel - is present, making the BBN bounds less stringent. 

In this section, we discuss the concrete examples of a relic stau or stop in the MSSM. 
We use the MICROMEGAS package to take into account all relevant annihilation 
and co-annihilation channels, but compare also with the results for Yab or Ynab for the 
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Figure 7: Maximal total yield YJ^^^ allowed by BBN as a function of the relic's lifetime 
Tx for the two cases Bh = 1 (red) and Bh = 10^^ (blue). The full lines are for a mass of 
mx = 100 GeV, while the dashed lines are for mx = 1 TeV. Note that for Bh = 10~^, 
the limit for tx > 10^ s comes from CBBN. 

case of one single gauge interaction. 

6.1 Relic stau 

Our results for an electrically charged relic can be applied, for instance, to the case of 
the supersymmetric partner of the r. We assume here that the relic stau is a right-chiral 
state, f/j, and that all other SUSY particles as well as the heavy Higgs bosons decouple. 

The dependence of the yield on the stau mass is shown in Fig. [HI For a 100 GeV 
fpi, we get Yf = 4.8 x 10^^^ at tree level from annihilation into photons (c.f. the dashed 
line). This is reduced by about 12% by the Sommerfeld enhancement (dashdotted line). 
In the full EW theory, the stau also annihilates into W~^W~ , ZZ and 7Z. In fact, for 
rrif = 100 GeV, the 77 channel contributes about 55%, •yZ about 25%, ZZ about 10% 
and WW about 5% to the total rate; the remaining 5% go into SM fermions. At higher 
stau masses, we have ~50% 77 and ~30% •yZ. Overall this gives a reduction of F by a 
factor of about 2 (sohd line), leading to Yf = 2.4 x 10~^^ at nif = 100 GeV. 

Staus can also annihilate into rr through t-channel neutralino exchange. We here 
consider only the bino contribution. Lowering the bino mass m ^ decreases the yield until 
bino-stau coannihilation takes over, increasing it again. We find a minimum yield at 
about rrij^ ~ (1.1 — 1.2)mfj, shown as dotted line in Fig. [8l It is roughly a factor 2 lower 
than the solid line, in agreement with [25j. Note also that the neutralino exchange leads 
to annihilation of same-sign stau pairs, f^f^ — > r^r^, so this process gets Sommerfeld- 
suppressed, and the total Sommerfeld effect almost cancels. 

The annihilation into W~^W~ and ZZ is considerably enhanced if the relic stau also 
has some fi component, fi = f/j sin 6 + fL cos 6 with cos 6^0. In this case also t-channel 
exchange of z/^ (for VT+VT") and f2 (for ZZ) has to be taken into account in addition 
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Figure 8: Total yield Y (black lines) for a relic as function of the stau mass. The dashed 
and dashdotted curves are tree-level and Sommerfeld- corrected results, respectively, from 
annihilation into photons (i.e. Yab)- The full line includes also annihilation into W and 
Z bosons, assuming all other sparticles decouple. Finally, the dotted line shows the case 
= lArrif-^. BBN bounds are shown in red: as full line for 0. 1-100 s lifetime and 
Bh = 0.65, and as dashed line for > 100 s lifetime and Bh = 10^^. Note that if the 
lifetime exceeds about 10^ s, the CBBN constraints become more important and quickly 
exclude number densities at the level of 10~^^-10~^^, see Fig. [71 

to the 4- vertex and s-channel 7/Z exchange, c.f. Appendix B. It turns out that these 
t-channel diagrams lead to a destructive interference: for given cos 6', smaller Uj. and f2 
masses lead to smaller cross sections. Since the stau and sneutrino masses and stau 
mixing angle are related to each other, one cannot simply maximise the cross section by 
choosing maximal mixing (cos^ = 0.7) and very heavy Ur and T2. However, for reasonable 
parameter choices, it is still possible to reduce the yield shown in Fig. M by up to about 
an order of magnitude. Alternatively, one could rely on resonant annihilation through 
s-channel Higgs exchange or on coannihilation with sparticles that are close in mass to 
bring below the BBN bounds. 

Barring these possibilities of largely enhanced cross sections, the stau lifetime and 
branching ratio into hadronic modes become key parameters to decide whether the sce- 
nario is allowed. First of all, let us discuss briefly the branching ratio into hadrons. We 
are considering here the decay — > r+LSP. The r decays into charged mesons 65% of 
the time, while the remaining times into leptons only. Charged mesons have a similar 
effect as nucleons during BBN only at short times < 100 seconds, because later they 
decay before interacting with nucleons and give rise only to electromagnetic showers [ID] . 
Therefore we will take Bh{t) ~ 0.65 for lifetimes up to 100 s, while it becomes much 
smaller for longer lifetimes, we will use BhIt) ~ 10^^ as reference value. This is in the 
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nif [GeV] [TeV] 

Figure 9: BBN-excluded regions for a gravitino LSP in the plane 777.3/2 vs. rrif. On the 
left a zoom on nif = 50-500 GeV for Yf = Yah (light grey ), = Yab/S (medium grey) 
and If = Yah/'^ (dark grey). On the right for = 0.1-10 TeV. Note that LEP excluded 
nif < 99.4 GeV for a charged particle stable within the detector [50] . 

central range computed recently for the stau decay into tau, gravitino and a qq pair, and 
we refer to that result for a more detailed analysis [H]. (A full computation including a 
more complete treatment of the hadronic decays of the tau for the case of a right-handed 
stau has been given in [52].) We have then to apply the BBN bounds discussed in the 
previous section corrected by these branching ratio factors, according to the time of decay. 

Regarding the stau lifetime, this depends strongly on the nature of the LSP. For the 
case of the axino LSP, the decay rate is given by 

s . ^ 1 9/ \f nir, \2/lO"GeV\V rn^A , . 

r(f„ ^ ™) = (25 (^) (^) j (1-4) (59) 

where is the axino mass, is the Bino mass, is the Peccei-Quinn scale, and ^ is 
a factor of order 1 taking into account some uncertainties in the loop computation [JU]- 
Therefore only the weakest BBN bound applies and actually disappears completely for 
large stau mass: in fact even for the conservative case 7/7^ = Wvtif and = 10^^ GeV, 
the lifetime becomes shorter than 0.1s for vnf < 590 GeV. We are here neglecting the 
case of a strong degeneracy between the stau and axino masses. We see therefore that for 
axino LSP a very light stau is a viable possibility and, depending on the supersymmetric 
spectrum, only the mass window between 125/250 — 590 GeV is possibly excluded by the 
BBN constraints, as can be seen from Fig. [HI In that region however probably a more 
proper computation of the stau hadronic branching ratio and its effect in the early stages 
of BBN is needed, as discussed in |12]. In fact comparing our exclusion region with theirs, 
we find that their constraints are much weaker for short lifetimes, due to an up-dated 
value of the Helium abundance and a larger systematic error, allowing all the stau region 
for an axino LSP. 
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For a gravitino LSP, the decay rate is given by 



5 /lOOGeV 

V "^3/2 



) 



2 




4 



r(fij ^ tG) = (5.9 X 10^ s) 



-1 



( 



100 GeV 



) 



(60) 



which typically gives longer stau lifetimes than the axino case. Figure shows the BBN- 
excluded region in the 7113/2 vs rrif plane. We consider a number density Yf equal to 
1/2 and 1/4 times Yab to account for the possible variation depending on m^. As can 
be seen, to avoid all bounds we need either a very light gravitino in the MeV range for 
rrif ~ 0(100) GeV, or a very heavy stau, e.g. rrif > 1.4 TeV (9 TeV) for 7/13/2 = 1 GeV 
(100 GeV), corresponding to a stau lifetime shorter than 0.1s. On the other hand, for 
ruf ~ 100-250 GeV and a lifetime longer than 100 s, Bh ~ 10~^ can bring the effective 
yield below the bound of mV ^ 5 x 10~^^ required by hadronic showers. Last but not 
least, note that the constraint from catalysed BBN becomes stronger than the hadronic 
ones for lifetimes longer than about 10^ s and excludes a light stau NLSP for gravitino 
masses above 10-100 GeV. 

6.2 Relic stop 

To discuss the case of a relic stop, we assume that only tpt is light while all other SUSY 
particles are heavy and decouple. Moreover, we assume that the light Higgs is SM-like 
with a mass of rrih = 115 GeV, and that the other Higgs bosons are also heavy and do 
not contribute to the stop annihilation. 

Results for the yield as a function of the stop mass are shown in Fig. [10 Let us first 
discuss the left plot, Fig. [TOlfa). which shows the yield at leading order (LO). Here the full 
line is the pure QCD result, Y^ab for SU{3), without Sommerfeld correction. As can be 
seen, tt* — >■ gg alone is efficient enough to avoid the BBN constraints up to stop masses 
of about 700 GeV. In the full theory, the stop can also annihilate into other particles, 
in particular into EW gauge and Higgs bosons. The yield for the QCD+EW case, still 
assuming heavy sparticles, is shown as the dashed line in Fig. [TOlfa). The dip at ~ 120 
GeV is due to the onset of t^t|j — > hh. Other important channels are annihilation into 
W^W~ and 7(7, contributing about 10% each to the total annihilation cross section for 
mi > 200 GeV. Annihilation into ZZ contributes about 5%. Annihilation into top quarks 
is suppressed by the heavy gluino mass, and also by rrit. However, if mi > 200 GeV 
and mg ~ 2mj, t/jt/? —>■ tt further reduces the yield by 10-20%. This is shown as the 
dash-dotted line in Fig. [TOT a). All in all, annihilation into gluons is, however, always the 
dominant channel, contributing at least 50%. We therefore take l^nafe/2 as a rough limit, 
which is shown as the dotted line in Fig. [TUlfa). Comparing with the BBN constraints we 
see that a relic t^ with a lifetime of 0. 1-100 s can be in agreement with BBN even for 
high masses of about 1 TeV. 

The impact of the Sommerfeld enhancement is illustrated in Fig. [TO](b) for the case 
tt* gg. As can be seen, taking the averaged Sommerfeld factor of C^^^-g^ = 11/42 in 
eq. reduces the LO yield by roughly a factor of 2, while a summed factor according 
to eq. (JMl) reduces the LO yield by roughly a factor of 3. These results are in qualitative 
agreement with those of [20], that considered the Sommerfeld correction in the neutralino- 
stop coannihilation region. Here note that for colour-singlet channels like, for instance, 
tt* —>■ W~^W~ a factor of C = 4/3 applies, hence leading to even larger enhancement. We 
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Figure 10: Results for for a relic tn as a function of the stop mass. In (a), tree- level 
results for different channels: the solid line comes from tt* gg only, the dashed line 
includes all channels into QCD+EW gauge and h bosons (case of decoupled sparticles 
and heavy Higgses), the dash-dotted line is the result for nig = ^m^^, and the dotted line 
the limit = Ynab/2. In (b), the effect of the Sommerfeld enhancement on the yield from 
tt* — s> gg: the full line shows the tree level result, the dashed line the result for a^^^^, i.e. 
applying an averaged Sommerfeld factor Cg^^^-^ = 11/42, and the dash-dotted line is for 
^sFsum^ i.e. applying a summed factor according to eq. ( IMI) . The BBN bound for 0.1-100 s 
lifetime is shown as thin red line in both plots. 



leave a detailed numerical analysis of the enhancement of the various stop annihilation 
channels for future work. Here we just note that the overall effect can be a reduction of 
the yield by an order of magnitude. 

Additional annihilation can take place after the QCD phase transition, when the stops 
are in a confined phase with the quarks. Since the lighter fermionic mesons are neutral 
and assuming that the annihilation process takes place without the formation of a bound 
state between the mesinos, the unitarity cross section is probably a good estimate of such 
annihilation and allows for heavier stops to be consistent with hadronic shower constraints. 
We see in fact from Fig. [H] that the unitary cross section with three degrees of freedom 
gives a yield well below all the BBN bounds (and below the range in Fig. [TOj) for stop 
masses up to 700 GeV. If also bound states between the mesinos can form efficiently, the 
BBN constraints disappear altogether [31], but note that we do not have to rely on the 
enhancement coming from such processes, which are very difficult to compute, for a wide 
range of parameter space. 

Let us briefly discuss the lifetime also for the stop case. For the case of an axino LSP, 
the stop decay rate is a larger than for the stau since it depends on the gluino mass and 
the QCD gauge coupling [3J: 

o X 1.9/ m \f m-a \VlO"GeV\V-, 

r««^«a) = (OxlO-secrt?(^)(^) j 

(61) 

where E,t is again a factor of order one taking into account the uncertainties in the loop 
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Figure 11: Lifetime of a relic t/? in the plane ^3/2 vs m^-. Recent results from the Tevatron 
exclude a metastable stop below 250 GeV [52j. 

computation [l9], in principle different than the one for the stau. Therefore, for the axino 
case, the BBN bound never applies if the decay into top is kinematically allowed, i.e. if 
mf > {ma + mt)'^. If the stop mass is smaller, the decay can proceed through a virtual top, 
for which we estimate a suppression of order (9(1/100) due to the 3-body phase space. 
This would still give a lifetime of order 0.1 sec, so the BBN constraints are completely 
avoided, as long as there is not a strong degeneracy in mass between LSP and NLSP or 
the factor C,t is exceptionally large. 

For a gravitino LSP, on the other hand, the same formula applies for stop as for stau, 
eq. (160|) with f — > t , because the gravitino couples only to mass. Note, however, that also 
in this decay the width gets phase-space suppressed if < mt + rriQ. For illustration, we 
show in Fig. [TT]the band of 0.1-100 s lifetime in the plane ^3/2 vs m^-. For lifetimes longer 
than 100 s, stops can still be in accord with BBN thanks to the additional annihilation 
during the QCD phase transition, if their annihilation reaches the unitarity one. We 
therefore conclude that cosmologically stops are an allowed NLSP in any mass range and 
in particular also for a heavy gravitino. Our results are in agreement with those for specific 
supersymmetric models with stop NLSP discussed in |51j. From the colliders side, note 
that the low mass region nii < 250 GeV has been recently excluded by the search for 
charged massive particles at the Tevatron |52j . 

7 Conclusions 

We have studied the number density of a charged relic by computing the annihilation 
cross section into gauge bosons, including the Sommerfeld enhancement. We have found 
that the Sommerfeld factor increases the thermally averaged annihilation cross section by 
20-50% and reduces the final yield even by a factor 2 or 3 for the SU (3) case. Moreover 
the result is very sensitive on how the higher orders are resummed. Nevertheless the 
number density surviving the annihilation is still large and BBN constraints are relevant 
for most relics. They can be avoided completely only for very large for particles in the 
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fundamental representation of SU{N) {N > 100 for mx < 10 TeV) or for cross sections 
nearly fulfilling the unitarity bound. For the cases of SM gauge groups, the allowed regions 
only correspond to very light relic masses, where the number density is low enough, or 
to sufficiently heavy relic masses so that the decay takes place in the first stages of BBN. 
The latter allowed region depends strongly on the relic decay channel, and, in case of 
a gravitino LSP with conserved R-parity, also on the gravitino mass. Let us mention 
here that if R-parity is just marginally broken, the NLSP can decay with shorter lifetime 
through R-parity violating channels and the BBN constraints can be easily evaded for 
any NLSP while keeping the gravitino LSP as Dark Matter |53j . 

More specifically, for the stau NLSP the light mass window has nearly completely been 
excluded by direct searches at LEP, even if the annihilation cross-section is maximal ~ 
4a{ff* 77), unless the gravitino is lighter than a few tens of GeV, while the large mass 
region is unfortunately out of reach at the LHC for gravitino masses m3/2 > lOOGeV. The 
detection of a quasi-stable stau at the LHC would then point to a scenario with relatively 
light gravitino mass, R-parity breaking or an axino LSP and could probably exclude the 
gravity mediated supersymmetry breaking scenario. In that case the determination of 
the stau lifetime and its decays will become crucial in distinguishing the different LSPs 
[71119]. 

The stop case is much less constrained thanks to the stronger annihilation cross- 
section, even if in this case the decay always produces mainly hadrons. We have practically 
no constraints if the LSP is an axino and even for a gravitino LSP, we can allow for 
relatively light stops up to approximately 700 GeV (1 TeV for lifetimes below 10^ s), if 
the annihilation cross section reaches the unitarity one after the QCD phase transition. 
The window between the present Tevatron bound around 250 GeV and 1 TeV should be 
surely completely covered by the LHC, the signature being a quasi-stable heavy fermionic 
meson. The detection of such a state would call for a non-minimal SUSY breaking sector 
with a coloured NLSP and a very weakly interacting LSP. In this case again only the 
analysis of the stop decays would allow to distinguish between the lightest states. 
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A Annihilation into massless SU{N) gauge bosons 



A.l Amplitudes for the annihilation 

We consider the case of one particle and antiparticle in the representation Tf and its 
conjugate, with momenta pi,p2 and mass m annihilating into two massless gauge bosons 
with group indices a, b, momenta P3,P4 and Lorentz indices /i, u respectively. 

The process has four different contributions, corresponding to the following four Feyn- 
man diagrams: 

t particle exchange in the t-channel described by the amplitude 



u particle exchange in the u-channel described by the amplitude 



(63) 



note that this contribution is identical to the t-channel under interchange of a ■(-^ 6, 



4 supersymmetric four-scalar coupling giving the amplitude: 

this contribution is symmetric in the exchange of a, b and therefore also /i, u; 



(64) 



s off-shell gauge boson in the s-channel decaying into two bosons via the non-abelian 
interaction: 



+(pi-P2)'^(2p3+P4)1; 



-[g^^-^{t-u)-{2p,+p3r{p,-p2y 



(65) 



this contribution is completely antisymmetric under the exchange of the gauge 
bosons group indices and therefore also under the exchange of their momenta and 
Lorentz indices. 

For convenience, we can then separate the amplitude into symmetric and antisymmet- 
ric part in colours a,b; then the interference between the two parts vanishes. Using 



rparpb f T'^ T^^\ -1- — 



(66) 



we have then 



^9% 



L J jt 



{2p^-p?y{2p2-p4Y , {2p^-P4r{2p2-PzY I y 



t — w? 



+ 



u — m" 



(67) 
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and 



rpa rph 



(2pi - P3)^(2P2 - ^ (2pi-p4)"(2p2-P3)^ 



t — w? 



u — 



(6^ 



In the Boltzmann equation, we have to insert the averaged cross- sect ion, so first we have 
to sum over all the final and initial states, i.e. sum over the gauge bosons polarisations 
and over all the group indices. 



A. 2 The matrix element 

The computation for the symmetric piece is straightforward: 



sym I 



9%\[t\T'}^^ 



2\2 



+ 



(u + m 



2\2 



+ 



Am 



2\2 



{t + m 



^ s/2 - 4m2 - 2(t - m^) ^ s/2 - Arn^ - 2{u - m^) 



t — 



+ 



2m^ 



(t — rrfiY t — m 



u — 

2im? ( 2m^' 
1 



(69) 
(70) 



In the antisymmetric part instead we have to take into account ghost subtraction and 
the total result is 



\A 



asym \ 



9n\ 



rpa rpb 



— + ^ - 4 



2\2 



1 is-Am^f {t-uf 16m2 
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So for the total matrix element we have 
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and the cross section is given in eq. ffTUj) . 



(71) 



(72) 



(73) 



31 



A. 3 Comparison with QCD result 

For the case of SU (3) we have 



E |{nT^} 

a,b,i,j 



28 

y 



and 
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a,b,c 



abc 



12. 



So after the sum over colours, we get 



\M\' 



14 [t-uf 48m2 2m^ 



+ 



2rn? 



28 



t — m? 



+ 12 



s 

2m2 



28 



{t - m?Y V 3 



28 

+ 12 

3 



+ 12 



^9t 



32 t-m^ [t-ni^Y 48m'^ 
— + 24 + 24^ + 



128 128 m2 

+ — ^ ^ + 



3 (t — m^)^ 3 t — 



1 + 



1 m 
iT 



This result coincides with the one given in the literature for the QCD case 
pare in general with [26] : 
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using again the symmetry in m ^ t and eliminating u. 
We have also that 

Co=Y: fabc = N{N' - 1) = 24 Ck 



N 



and for a single RH stop, we must use 2nf = 1. Then we get 



iMigg^tRtji) 
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which coincide with our result above eq. (177|) . 
Now integrate over t and obtain 



a m, s) 




(82) 



which coincides with [2S] allowing for the exchange of initial and final state (s — 4m^ s 
in the denominator) and the initial state averaging, i.e. a factor of 1/64 for the two gluons 
initial state. 



B Annihilation into SU{2)l gauge bosons 

Another important channel of annihilation for light stops or staus is into EW gauge bosons. 
Let us consider first the pure SU{2)l case, neglecting the gauge boson masses, but with 
a split SU{2) multiplet. We consider here the case of one left-handed sparticle and one 
left-handed antisparticle of momenta pi,P2, mass mi and SU (2) index 1, annihilating into 
2 gauge bosons of SU{2)l index i, j, momenta P3,P4 and Lorentz indices /i, v respectively. 
Then we can directly use the result for SU (N), only taking into account that —>■ 
with (T* denoting the Pauli matrices, and that in this case we have an initial state made 
of the upper components of the SU{2)l doublet, while the lower component is exchanged 
in the t- and u-channel and can have a different mass m2. 

We have then for the two amplitudes, symmetric and antisymmetric in the group and 
Lorentz indices, 
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To compute the annihilation cross section, we have to sum over all the final states and 
initial states; this means that we have to sum over the W polarisations and over the 
SU{2)l indices but in this case the initial state group indices are fixed. 
The symmetric piece gives 
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In the antisymmetric part instead gives 
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B.l SU{2)i sum and total matrix element 

In this case the sum over the indices i,j is simple. We have that 



^ 4 



J = Eii2^iAir = E^i = 2+i 



(89) 



where we have considered the annihilation into W^''^ separately from that into W3. In 
fact the intermediate particle has a different mass in the two cases. 
On the other hand the antisymmetric product gives 



(90) 



since in this case only can be exchanged in the s-channel for W^'"^ in the final state. 
Then the matrix element for annihilation into W^''^ gauge bosons is given by 
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while the annihilation into has only the abelian contribution with the presence of a 
single mass 7771 
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The cross section for the first case is then 
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while the annihilation into is identical to the abelian one in eq. ^ for ex = 1/2. 



(92) 



B.2 Annihilation into physical W^W 

Let us now consider the case of a broken SU{2)l symmetry like the Standard Model 
and massive gauge bosons which mix to give the physical ,W~ , Z,'y. At the same 
time let us consider a general initial state given by the light stau mass eigenstate fi = 
tl cos Of + fn sin 6f and its antiparticle. In this case the intermediate particle exchanged 
in the t- and u-channel can be only a left-handed sneutrino and therefore we can neglect 
the mixing for the intermediate state. 

Then the annihilation into W~^W~ is given by the following channels: 



t sneutrino exchange in the t-channel described by the amplitude 
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(93) 



u NO u-channel since and W are different particles ! 
4 super symmetric four-scalar coupling giving the amplitude: 
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this contribution is symmetric in the exchange of /i, z/; 
s off-shell Z/'-y in the s-channel decaying into two WW via the non-abelian interaction: 

1 



■9l 2n (. 4sin2 0vv' 
I— COS Of\\ — 



3 cos2 Of s- Ml 



[g^''{t-u) 



(95) 



-{2PA + P^)^{pi - P2y + {Pi - P2)^{2p^ + PaY 
+ie'\- [g^''{t-u) - {2pA+P3)'{pi-P2)'' 

+ (P1-P2)^(2P3+P4)1 

.gl 2 / AsmHwMl\ 
cos Of \ 1 ' 



2 ' \^ 3 cos^ Of s J 

g'"'{t -u)- {2pA + P3)^{pi - P2Y + {pi - P2)^{2p3 + Pa)' 

s-Ml 



(96) 
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this contribution is completely antisymmetric under the exchange of the W momenta 
and Lorentz indices. Note that the photon contribution is proportional to = 
gi| sin^ 9w and cancels exactly with the second term due to the Z-boson in the 
case of equal mass. In that limit in fact the U{1)y factor decouples and does not 
participate in the non-abelian interaction. 

s-H off-shell h/H in the s-channel decaying into two WW via the non-abelian interaction; 
in this case we have to consider both neutral Higgses: 



sH 



2 ^ s 



Chs 
s-Mjj 



+ 



Ml 



(97) 



where Cu/h is coming from the product of the coupling of the staus to the Higgses 
and of the Higgses to the WW pair. These constants depend on the whole SUSY 
breaking parameters. For the staus these couphngs are probably neghgible. We 
have in fact 



C 



H/h 



4 

- sin^ Owiy — tan^ Of] 



(1 +tan2/?)cos4 Ow 
+^ «,2cos2^^(l + tan2/?) + 



, fi ^{Z\H/h + Z2H/h tan/3) 
g2 cos^ 9wMw\i + tan^ (3 

X [Z2H/hAr + Z^H/hiK + A^*^r) + h.C.) , 



(98) 



where Z is the matrix which diagonaliscs the Higgs mass matrix, is the tau 
Yukawa coupling, A^,A'^ are the SUSY breaking trihnear terms and /j, the Higgs 
supersymmctric mass parameter. This contribution is suppressed by M'^/s for 
large s. We can include it easily into the 4- vertex contribution by substituting 



w 



Ml 



Kh{s) . 



(99) 



Now we can write the t-channel as the sum of a symmetric and antisymmetric part, 
adding and subtracting a fictitious u-channel, as 



(2pi-p3)^(2p2-p4)" 

t — ml 



1 
2 

1 

+ 2 



+ 



(2pi -p^Y{2p2-p^ Y , (2pi-p4)"(2p2-P3) 
t — ml 

{2pi - P3Y{'^P2 - PaY 



u — m~, 

(2P1-P4)"(2P2-P3) 



t 



m~ 



u 



ml 



(100) 



so that we can make contact with the previous computation and find for the symmetric 
and antisymmetric amplitudes respectively: 



q2 

-Fi^cos^ 9f 



{2p^-p3Y{2p2-Piy 



(2pi - P4)"(2P2 



t — ml 

PsY 



u 



ml 



+ 2g^^{l + KH{s)) 



:ioi) 
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and 



■92 2 n 
I— COS U:f 



t-mi 



(102) 



{'2pi-P4y{2p2-P3)^ 



+ 2{1-Gz{s))x 



X 



u — ml 

gi'^t -u)- {2p^+p^Y{pi -P2T + (pi -P2Y{2pz+PaY 



4 sin^ Bw 
3 cos^ 9f s 



s-Ml 



vanishes in the hmit of zero Z mass. This coincides with 



where Gz{s) 

the previous result for Kh, Gz, Mz = 0, a part for a sign, which just corresponds in 
exchanging i j. 



B.3 Polarisation sum 

The sum over the W polarisation in this case is given by the polarisation tensor 



_ _|_ 



a u 
P3P3 



(103) 



where pa is the gauge boson momentum. 
We have then for the matrix element 



\f3A 



/if I 



HV I 



Ml 



Ml 



+ 



Ml 



(104) 



in this case neither amplitude vanishes when contracted with the gauge boson's momen- 
tum. Note that the second and third contributions are related again by the symmetry 
P3 <-> P4; <-> // and are equal since the final state has two particle with the same mass. 



B.4 Symmetric part 

We must compute the four contributions, and we have then 

^Icos^^, r , {ml + mj-Ml/2y 



A* _ 



1 ■ 



+ 2Kh{s){1 + Kh{s)) (105) 



2 f ' {t-mlY 
{s-Aml + Mlf 



s + 2(m? - m? - Ml) 
cos^ Of 



+ Kh{s){s - Ami - + 2M, 



wJ 



i + M±!!i^W^2i..(.)(i + i.,(.)) 



{t - mlY 



1 ( 2^0 2 .,.2 1 {2ml + 2ml-^Mlf 
■^——2 I + 3m^ - 2M^ 



2 s + 2(m? - m? - Ml) 



Kh(s) 



{s - Ami - + 2M, 



w) 



(106) 



which in the limit of vanishing M^ and = mf coincides with our old result. 
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The other pieces give instead 



IJ.!/ I 



+ 



\plA 



jXV I 



cos^ 6'f 



2\2 



2ml + 2m| 



t-mls + 2{ml - m? - M^) 



ml 



ml + s/2" 



ml 



(107) 



2Kl{s) 



and the last part: 



fif| cos^ Of 



m. 



2\2 



2\2 



^^ml - m| 
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o 

m 



'l-l 



1 + 



1 {ml - m|) 

2 (t-m?)2 



2 2 

mr, — ml 



As + 2{ml-mf + M^) 



s 
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- 2 



(108) 



Both these contributions vanish in the hmit of equal stau and sneutrino masses and zero 
gauge boson mass as they should. 

So summing all together the result is 
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+ 



2ml + 2ml - Mil, 



s + 2{ml 



2 

m~ 



Mir) 



m 



2\2 



Ml, 

2 2 

mi — m| 
t — ml 



1-1- 



2 (t-m?)2 

2 2 
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2 2 
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2MI,' 
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(109) 



Note that the in the limit of large s, sK]-[{s) remains finite and therefore there is no 
problem with unitarity. 
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B.5 Antisymmetric part 

The antisymmetric piece is more involved. We have 



A* A^"" 



cos^ 9f 



t — m~ 



2 + 2 



[ml + m| - M^/2)- 
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4(m? + m?) - 2M^ + 



[s - Ami + 



2 \2 



+ {l-Gz{s)) 

+ (1-Gz(5)) 
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M|)(t 



m- 



4(s + m? — ml){s — 4m? 
(s-M|)(t-m?) 



which in the hmit of vanishing My[r, Mz and m^ — mf coincides with our old result 
The other pieces give instead 
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'-fiv I 
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(111) 
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in the hmit of vanishing m? — m?, M^, masses keeping Mz/Mw — > 1 we have 
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Ml 



§2 COS^ Of it — u)' 
4 27~ 



(112) 



as expected from the QCD result. 
The last part gives instead 
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g-i cos 
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IQMl 
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(113) 



+ 2- 



(t-ml) \{t-ml) s + 2{ml-ml-Ml)J_ 
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Note that this contribution does not vanish in the hmit of equal stop and sbottom masses 
and massless gauge bosons. In fact keeping Mz/Mw — > 1, we have 



HI/ 1 



8 



2s 



(114) 



which gives the annihilation into the Goldstone part of the Higgs field. 

We can now put all together to give 
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Note that to reduce these expressions in terms of only the t variable, we have used the 
simple decompositions, i.e. from s + t + u — 2m? + 2Ml one obtains 
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B.6 Results for the cross section 

We can integrate the matrix element to obtain the cross section in the two cases: 
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